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A fully algebraic approach to reconstructing one-dimensional reflectionless potentials 
is described. A simple and easily applicable general formula is derived, using the 
methods of the theory of determinants. In particular, useful properties of special 
determinants - the alternants - have been exploited. The main formula takes an 
especially simple form if one aims to reconstruct a symmetric reflectionless potential. 
Several examples are presented to illustrate the efficiency of the method. 
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I. INTRODUCTION 


Formulation and solution of inverse problems is an increasingly important field of scientific 
research. However, compared with a well-posed (in Hadamard’s sense) direct or forward 
problem, the corresponding inverse problem is much more difficult and, as a rule, ill-posed. 
Inverse scattering problem can be considered an exception to this rule. Namely, in the 
simplest one-dimensional case, the inverse scattering theory provides strict mathematical 
criteria for the existence, uniqueness and stability of the solution. It means that in this 
particular case the inverse problem is well-posed as well. 

The related forward problem is the solution of the simplest time-independent Schrodinger 
equation 



( 1 ) 


for a given potential U(a;), and subjected to appropriate physical boundary conditions. Eq. 
dU) can be easily solved numerically and thus, in principle, all spectral characteristics of the 
potential V ( x ) can be accurately ascertained. 

The inverse problem is to determine the unknown potential starting from the known 
spectral characteristics. This is a serious task even in this simple case for the following 
reasons: 

1. It is not obvious what kind of input information is actually needed to solve the problem 
uniquely. 

2. There must be a theoretical basis (a fundamental equation) which enables to solve the 
problem. 

3. Apart from the theoretical difficulties, another important question arises: how to 
obtain the necessary input data? 

4. Even if the mentioned principal barriers could be overcome, the computational- 
technical solution of the problem is not at all trivial. 

Issues 1 and 2 have been successfully resolved in early 1950s for the class of potentials 
on the half line: x 6 [0, oo). The necessary and sufficient conditions for the unique solution 
of the inverse problem have been formulated in a series of outstanding theoretical works 
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by Marchenko^, Gelfand, Levitan^, KreinA’^ and others (see, e.g.,-, section III.7 for a brief 
overview). In addition, three different methods to solve the problem have been worked out, 
based on the integral equations by Gelfand-Levitan^, Marchenko^ and KreinA 

Solution of the inverse scattering problem on the full line (—oo < x < oo) is a more 
challenging problem, which was first addressed by Kay-, Moses^, and few years later in a 
series of papers by Faddeev^— . These fundamental studies provide full description of the 
solution procedure, while the correct necessary and sufficient criteria for the uniqueness of 
the solution were given by Marchenko (see^, section III.5). These criteria apply to the 
following pair of integral equations&AEdih 


K 1 (x,y) + Ai(x + y)+ / A^z + y)K 1 (x, z)dz = 0, 


( 2 ) 


K- 2 (x,y) + A 2 (x + y) + / A 2 (z + y)K 2 (x, z)dz = 0, 


( 3 ) 


which are formally very similar to the Marchenko equation on the half line (see^, p. 218). On 
the other hand, their operator-theoretical content is closer to the Gelfand-Levitan approach. 
Therefore, as a kind of compromise, Eqs. (I2l) - (l3l) are often called Gelfand-Levitan-Marchenko 
(GLM) equations. The kernel A 1 in Eq. (J2J) (where y > x) is completely specified by the 
so-called right scattering data, while the kernel A 2 in Eq. (J3]) (y < x) is determined by the 
left scattering data. These two sets of input data are equivalent: the left data are uniquely 
determined by the right ones and vice versa. In the following analysis we will rely on Eq. 
(|2l) . If one is able to solve this integral equation then the potential is given by 

f dKi(x, x ) 


V(x) = -2 C- 


dx 


( 4 ) 


Now, let us briefly discuss Point 3 in the above list. Unfortunately, the criteria for the 
uniqueness of the solution are so strict that it is nearly impossible to get all the necessary 
input data experimentally. However, if one sets an additional constraint that the resulting 
potential must be reflectionless then the inverse scattering problem can be solved much 
more easily. Moreover, a symmetric reflectionless potential is uniquely determined if its 
full spectrum of bound states is known. In the following analysis it is assumed that the 
potential V(x) that corresponds to Eq. (J3J) is reflectionless by definition. 

As the general principles of building confining reflectionless potentials are long known--’- 4 , 
a natural question arises: is there any need to revisit the topic? A motivation comes from 
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Point 4 stated above: if the total number of bound states is large then the technical side 
of the procedure becomes important. This in turn motivates development of more efficient 
algorithms. In this paper, a new analytic algorithm is derived, which enables to easily 
calculate a reflectionless potential with an arbitrary number ( N ) of bound states. 

The paper is organized as follows. In Sec. II, the general principles are briefly described, 
which form the overall basis for the approach. Secs. Ill and IV make an excursion to the 
theory of determinants, the benefits of which are described and illustrated in Sec. V. Finally, 
Sec. VI concludes the work. 


II. RECONSTRUCTION OF REFLECTIONLESSS POTENTIALS: 
UNIVERSAL RECIPE 


Suppose we are given 2 N parameters for an unknown reflectionlcss potential V(x): the po¬ 
sitions of N discrete energy levels E n = —Ck^ and N norming constants C n (n = 1, 2,..., N ) 
for the Jost solution \l/ 1 (i«; n , x) of Eq. ©■ so that —» exp(— n n x) as x —* +oo, 

and 

OO 

I 'ff 2 n (x)dx = 1, T n = C r A> 1 (iK n ,x). (5) 

— OO 

Then it can be shownH that 


1 det 

= - yww 

Here A is a symmetric matrix with the following elements: 

, _ r | 

-ft-mn Omn i ? 

K m + K n 

A (n) 

is obtained from A by replacing the nth column with its derivative, and 


( 6 ) 

(7) 


A n (x) = C n exp(-K n x). (8) 

A simple formula can also be obtained for the potential^— 

V(x) = -2cA {In [ det (A)]}, (9) 

which is uniquely determined by parameters K n and C n . 

Note that the potential remains unchanged if det (A) is multiplied by a function exp ( ax + (5), 
where a and /3 are arbitrary constants. Consequently, we can multiply, for example, any 
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row of the initial matrix by a function exp [n n (x — x n )] (n = 1, 2, N), where x n are new 
parameters, equivalent to norming constants: 

C 2 

exp (2tt n x n ) = ——, (10) 

and thus, according to dehnition ([HD, 

^ = exp [~2n n (x - x n )]. (11) 

Z K; U 

As a result, we get another matrix which contains full information for reconstructing the 
potential according to Eq. (0: 


V(x) 


An — Bn + Cn , 

~ 2C i^{' n [ det ^ jv ) 


( 12 ) 


Here 


/ gKl(x—Xl) Q 



V 


o 

o 

o 


e K2(x-X2) 


0 N 
0 


(13) 


and 


/ g —Nl(®—Xl) 



2\/ k i k 2 g —ki fa;—xi) 
K1+K2 


2 v /fllK2 K2.[x — X2.) 

K1+K2 

e ~K2(x—X2) 


2^/KlK N ^ — K n (x — XKl) \ 

Kl+KJV 

2 ^k, 2K N ^_ Kri ( x _ XN ) 
K 2 +KJV 


(14) 


2^/k ikjv _ Kl (x—xi) 2 s /k, 2 Km K9 ( x —x 9 ) .-sjv(i-ijf) , 

\ K1+KJV C K2+KJV ... C / 

Remark: A subscript was added to denote the number of the bound states (and the 

rank of the matrix). 

From now on, the determinants having the structure det(A) • exp (ax + (3) with the ele¬ 
ments A mn defined in Eq. (!T|h will be sometimes called r-functions (as is common in soliton 

theory). In the next two sections it will be shown that such determinants can be easily 

calculated even for an arbitrarily large N. 


III. GENERAL FORMULA FOR THE r-FUNCTIONS 

To evaluate a non-trivial determinant, one can use the Laplace expansion (see^, p. 487) 
in terms of the fixed row (or column) indices. For example, choosing a set of indices mi, 
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777 - 2 , ■ ■■,iri k for an arbitrary N x IV-matrix An, so that 1 < m-i < m 2 < ... < N, we get 

det (An) = det A (mim 2 ...m k \n\ri 2 -..n k ) X (15) 

1 <ni <7i2 <... < AT 

det A (mi?77-2...mfc| n\n 2 --.n k ). 

Here H (mim 2 ...mfe| n\n 2 --.n k ) is a k x fc-submatrix of An that lies on the intersection of 
rows m 1 , m 2 , ■ ■■,m k and columns n \, 77 , 2 , ...,rik, while 

det H (m 1 m 2 ...mfe| n\ri 2 ---n k ) = 

( mim2 ... mfc | nin 2 ..m fc ), 

and M (m 1 m 2 ...mfc| 77i7i2..mjt) is a minor obtained from det (An) by deleting rows m-i, 
m 2 , ■ ■■,m k and columns 711 , 77 2 , ■ ■■,n k . 

Consequently, applying Eq. (TT5|) to the r-function defined by Eqs. (IT2]) - (fbi)) . 

N 

t n = det = a 0 exp(ao) + y a*exp(cq) + (16) 

Z=1 

y a ii exp(a ij )+ ^ a ijk exp(a ijk )+ 

l<i<j<N 1<Z<J<AJ<AT 

... + Gq23...Arexp(a 12 3...7v), 

where the coefficients a-o, a,, a^, dij k ,... as well as the corresponding arguments of the 
exponents can be easily fixed with the help of Eqs. (TT2]h(fT4]h Indeed, 

N N 

*) = !>(*- £j), a, = y (-l) 5i, /q(a; - xj), (17) 

1=1 1=1 

N 

Oij = y (~l) Sim+Sjm K m ( x - X m ), 

771=1 

N 

a ijk = y (-1 ) Sin+S ^ +Sk "K n (x - X n ), 

71=1 

N 

-01123...N = ^(-1 ) Sln+S2n+ - SNn K n (x - X n ) = —« 0 , 

71=1 

Su being the Kronecker’s symbol. 

To further simplify Eq. fllbp . let us group the terms into pairs, so that the arguments of 
the corresponding exponents differ only by sign. For example, the first pair is formed of the 
terms with coefficients 

a 0 = 1 
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and 


«123...iV = A 


2 _ 
0 — 


2 2 x /k 1 k 2 

Ki+K 2 

2y/~Klhi2 

«l+«2 


/■viH-Avjv 

2y/K2KN 

^2+^iV 


27^17 2 yj K2&N ^ 

^1+^iV ^2+^iV 

Note that 

det (Ejv] = a 0 exp(a 0 ), det (pN^ = Oi 23 ..jv exp(-a 0 ). 

Here we defined a new coefficient Ho, whose subscript ” 0 ” emphasizes that the expression 
—«o = —K\{x — X\) — hi2(x — X2) — ... — hi^(x — Xn) contains no terms (0 terms) with plus 
sign. The same logic can be applied to all terms of Eq. (fT6j) . For example, the appropriate 
partner for the term aiexp(ctq) is 0234 ..,jv exp(—an), where 


a 1 = 1, 


^ 234 ...N = A-l — 

-y 2 ^/k 2 K 3 2y/K2 K-N 

K 2 +K 3 K 2 + KJV 

2 v /K2^3 2 x /K3K N 

K2+K3 K3+KN 


2y/K2HN 2 X /K3KJV 
K 2 +KJV K 3 + KJV 


Analogously, we can form a pair from ai 2 exp(ai 2 ) and CI345...7Vexp(—a 12 ), where 


012 


1 

2 x /kiK2 

K 1 +K 2 


2 x /kiK2 

K1+K2 

1 


> 


2 x /K3Kjv 
K 3+KiV 
2y K 4 KN 
K,£-\-K,n 

2y / it3 Kjv 2 x /fC4Kjv 
«3+ k JV K 4 +KJV 

The principle is simple: for any a ni „ 2 ... ni exp(a nin2 ... ri J the partner is a ni+ini+2 ... njv exp(-a nin2 ... n 
In addition, as we will see below, it is convenient to define a relevant coefficient 


0345 ...iv — 


-j 2yK3M 

K 3 +K 4 


2y>i3K4 | 

K 3 +K 4 


d 2 

- rL n\ri2...ni 


— 0>n\n2...rii ' 0,n i+ 1 n i+ 2...n N , 


( 18 ) 
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where the indices point at the terms with plus sign on the right side of the expression 

^n\) "h ••• (19) 

+ K ni (x -X ni )~ K ni+1 (x - X ni+1 ) 

Looking at the structure of the matrix Am, it is obvious that all these plus sign terms 
can only originate from the expansion of det (pN^j and they correspond to the product 
(B]sr) nini ■ (. B]\r) n2 n 2 • ••• • (-Bjv) niTlj . The terms with minus sign on the right side of Eq. (jT9l) 
are related to the expansion of det (pN^ , without any contribution from det ■ 

On the basis of the above arguments, the following conclusions can be made: 


• All terms on the right side of Eq. (fT6l) can be grouped into pairs. There is only one 
term, 0123 ..,jv exp(—a 0 ) = det (pN^j (with partner det (pN^j = «oexp(a 0 )) which is 
entirely formed of the elements of the matrix Cm- Any other term (both partners) 
contains some diagonal elements of the matrix Bm as well. 


• Any term a nin2 ... ni exp(a nin2 _ n .) (i = 1, 2,...) can be obtained by replacing all elements 
of the rows and columns n l+ i, n l+2l ---TIm of the matrix Cm with corresponding elements 
of the matrix Bm (mostly with zeros). As a result, one gets a modified matrix Cm, while 
anin 2 ...ni exp(a niri2 ... r i i ) = det {pPj ■ The Laplace expansion (1T51) of this determinant 
for the fixed rows 7ij+i,7ii +2 , ...Un contains only one term! 

• For any term a nin2 ___ ni exp(cc niri2 _ n .) of the expansion (IT6|) there is a partner 




exp(-a nin2 ... n .) = det (C N ) , 


where Cn can be obtained by replacing all elements of the rows and columns 711 , n 2 , ■■■rii 
of the matrix Cm with corresponding elements of the matrix Bm- The Laplace expan¬ 
sion of det (Cm) for the rows 77 , 1 , 77 - 2 , ---rii also contains only one term: 


det (C^) = exp(-a 12 ...i)x 


2 \/ Kn i+l Kn i+2 

Kn i+l+ K "i+2 


‘ 2 V Kn i+l Kn i+‘2 

K n i+1 +Kn i+2 

1 


S+i+ k "jv 
V: n i+2 Kn N 

Kn i+2+ Kn N 


2 V K ’ i i+l Kn iV 2 \Z Kn j + 2 Kn iV 
Kn i+ l + K n N K ™i+2 + Krl N 


( 20 ) 
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It means, for example, that 


A 


2 

12 ...i 


— a l2...i ' Oi+l,i+2,...JV — 


2 2 ^/ 

K1+K2 

2^/K\K,2 | 

K1+K2 


2 y! K\Ki 2y/K2 Ki 
Kl+Ki K2 + Ki 


j 2^/Ki + lKi + 2 

Ki + l+Ki+2 
2\/'U + 1'm + 2 -j 

fti+l+ K t+2 


K2 + /^i 


X 


1 

2y/ K i+i K N 

Ki_)_l+Kjv 

-\f K - n i+-2 Kn N 

K "i+2 + K "iV 


( 21 ) 


2 V /Ki+1KJV 2 x / K n i + 2 K, n jv 

^i+1 “^^IV 


• One should avoid re-use of the terms: an already existing pair must not be included 
again! It means that the members of the modified expansion (fl6l) are identified by no 
more than [IV/2] indices (square brackets denote integer part of N/ 2). It is convenient 
to group the members on the basis of the number of indices, so that there will be 
0,1,2,..., [IV/2] different indices. 


As the final result of the above analysis, we get the following general formula: 

N 

t n = det ^ A [exp(c^ + $) + exp (-a* - /%)] + 

i=0 

••• + ^ [exp(Q!j 1 ...i [JV/2] + Pi 1 ...i [N/2] ) + 

l<«i<...<i[jv/ 2 ]<A r 

N 

exp(-a il .. i[JV/2] - Ai...^,)] = 2 A cosh (a, + &) + 

i=0 

2 Ai * 2 cosh(a'j 1 j 2 + /^jj) +...+ 

^ ] 2 ^-U*2...j[AT/2] cosh( 01 ^ 1 ^ 2 ...i[iv/2] A fiili2---i[N/2] ), 


( 22 ) 


where 


exp (A 


nin2--.rii) 


^n\n2...nj 
a ni +1 n i+ 2—nN 


(23) 


It is easy to be convinced that the expansion (122|) contains 2 N 1 members in total (apart 
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from inessential factor 2). Indeed, according to Newton’s binomial theorem 



where we took into consideration that only half of this formal series is actually needed. 


IV. ALTERNANTS OF r-FUNCTIONS 


We have shown that not only det itself but also the coefficients A ili21 ..., A ili2 _ i[N/2] 

in the expansion (1221) are r-functions. Consequently, the solution of the inverse scattering 
problem has been reduced to evaluating a number of determinants 


D(n i, k 2 , K n ) 


1 

2 s /k 1 K 2 

K1+K2 


2y^iK2 

K1+K2 

1 


2\/ N1 Kn 

2\/ N 2^n 
1^2 ~\-Kn 


(24) 


2-y/ 2^/ 

Kl+Kn N2+Rn 


fixed by the parameters Ki < K 2 < ... < K n , n being an appropriate natural number. We 


are now going to derive a simple formula for calculating such r-functions. First, let us set a 


one-to-one correspondence between each row of the determinant and a fixed parameter 


qt = (* = 1) 2 , -, n ) ■ 


(25) 


For example, the modified elements of the first row of Eq. ()24l) will be 


1 = 


Zqiy/Kl 2 ^/K x K 2 _ ^Qly/^2 ^y/^1 K n _ 2 q ls /K^ 


q{ + Ki Ki + K 2 qt + K 2 «1 + K n qt + Kn 

The usefulness of this trick soon becomes evident, although there seems to be only a formal 


change: 


D{K-L,K 2 ,..,K n ) = 


2 n ■ KiK 2 ...K n 


2 q\y/K\ 2qi^/K2 2gi x /«y 


qf+Ki 

2q2y/Kl 

?2+ k 1 

<?1+K2 

2q2y/K2 

<?2+ k 2 

9l + Kn 

2q2y/ 

q‘2~\-Kn 

2 Qny/fcl 

(j2+Ki 

2 Qn\/ K 2 
9n + K2 

2qnV^ 

ql+Kn 

1 

1 

1 


q'(+K i 
l 

91+K2 

1 

IJl+Kn 

1 


<?2+ k 1 

<?2+ k 2 



1 

1 

i 


9n + K l 

9n+«2 




(26) 


10 














































Here we separated a factor 2 q i: from each row and y HFi from each column, i.e., 2 qi^fFi = 2/tj 
(f = 1, 2, n) from each such pair. 

Next, let us transform Eq. (127)j) into a polynomial, multiplying each row by 


II(^) (i = 

j 


The result is 


hi 2 , ^n) — __ 2 ^ ^11) (27) 

II (Ki + Kj) 

l<i<j<n 

where a new determinant 


D 


n 


n (?? + n (e? + kj) ... n («? + «*■) 

jA2 


n (?2 + K i) n (?2 + ••• n + «*) 

jA2 j^n 


(28) 


nvi+ '.j i n («s+ K i ) n (<&+ 

jAl j¥2 


was introduced. As can be seen, the factor 2 n • aci «2 ■ ■■Av n was canceled out from Eq. (f27l) . 

We can see that the elements of the columns of D n correspond to different values of the 
same function, while any row is characterized by a single fixed parameter. Indeed, Eq. (I28|) 
can be expressed as 


*i(?i) F 2 (Qi) ••• F n (qi) 
Fi(q 2 ) F 2 (q 2 ) ■■■ F n (q 2 ) 


Fi{q n ) F 2 (q n ) ... F n (q n ) 


where 


F\(x) = ( x 2 + k 2 ) ( x 2 + k 3 ) ... (a^ 2 + «n) , (30) 

F 2 {x) = (x 2 + k i) (x 2 + hi 3 ) ... (x 2 + k„) , 

F 3 (x) = (x 2 + Ki) (x 2 + k 2 ) (x 2 + k 4 ) ... (x 2 + /c n ) , 

F n (x) = (x 2 + Ki) (x 2 + /C 2 ) (x 2 + K 3 ) ... (X 2 + K n _i) . 

A determinant that has a structure of Eq. (|29|) is called alternant (see^, p. 161). The best 
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known alternant is Vandermonde’s determinant (for the same set of variables) 


1 qi q\ ... qr 1 

1 Q2 ql qr 1 

V n = . 

i q n -1 ql -1 ••• ql~ i 1 

1 g n q 2 n ... g"" 1 

which can be easily evaluated (see^, p. 16 ): 


( 31 ) 


V n (qi,q 2 ,-,q n ) = {Qj ~ Qi) ■ (32) 

An important point is that the factor V n (qi, q - 2 ,..., q n ) can be separated from any nth-order 
alternant. Indeed, the argument q n may only appear in the nth row of Eq. (1291) : if we put 
it into any other row then the determinant would be identically zero. It means that D n has 

n— 1 

a factor n (Qn - qi) . Analogous reasoning applied to g n _i shows that D n also has a factor 
1=1 

n —2 

]1 {qn-i-Qi), etc. Putting it all together, we conclude that an nth-order alternate always 
1=1 

has a factor V n (qi, q 2 ,..., q n )- 

To continue the analysis, let us recall some useful properties of the elementary symmetric 
functions: 


(T 0 = 1, 


(33) 


&i = qi + q 2 + ••• + q n , 


°k= ^2 qhq% 2 ■ ■ ■ qi k , 

crn = qiq2 ■ ■ ■ q n - 

Here (Jk (k > 2) is a sum of all possible products of exactly k variables arranged in ascend¬ 
ing order of their indices. According to the Fundamental Theorem for symmetric polyno¬ 
mials (see^, p. 312), any such polynomial can be uniquely expressed as a polynomial in 
oq, <72, a n- This in turn is a basis for the following important theorem: 

Theorem 1: Let \A n \ be an nth order alternant generated by the functions 

Fj(x) = a 0 j + aij ■ x + a 2 j ■ x 2 + ... + a n j ■ x n , (34) 

3 = 2 ,ri, 
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where the parameters a t] do not depend on x, and define 


Then 


S k 


\An\ 

V n {.qi,q 2 ,-,q n ) 


-l) fc 

Vk{q 1 

, 92 ," 

1 9 n) • 

Ooi 

an 

a 2 i 

... a n i 

«02 

®i 2 

a 2 2 

... a n2 

O'On 

ain 

®2n 

• • • a nn 

S n 

Fn —1 

Fn —2 

... ^0 


Proof: Let us introduce an auxiliary polynomial 


( 35 ) 


(36) 


F(x) = (x- qi)(x - q 2 )...(x - q n ) 


(37) 


and form an (n+ l)th order Vandermonde’s determinant, adding a new (arbitrary) variable 
q n+ \, so that 


^ / n+l(?l) 92; ■■■) Qni 9n+l) 


1 

1 . 

.. 1 

1 

9i 

92 • 

•• 9n 

9n+l 

gr 1 

92 _1 • 

.. g n_1 

•• SZn 

9„;! 


92 ■ 

.. q n 

^n 

9n+l 


(compared with Eq. (T3TT) . the row and column indices are interchanged). Multiplying 
A n+ i by V n+i and using Eqs. (J33D, (J35]) , (J37|), we get 


Fi(qi) ^ 1 ( 92 ) 
F 2 (qi) F 2 (q 2 ) 


A n+ i X Vn+l — 


^( 91 ) F n (q 2 ) 
F(qi) F(q 2 ) 


Fi(q n ) Fi(q n+ i) 
F 2 (q n ) F 2 (q n +i) 

Fn^Qn) F n {c[ n -\- 1 ) 
^(9n) F(q n+1 ) 


^i( 9 i) ^1(92) ••• ^(g n ) 7 fi(g n+1 ) 

^2(91) F 2 (q 2 ) ... F 2 (g n ) ^ 2 (g n+ i) 

Fn( 9 l) F n (q 2 ) ... En( 9 n) F n (^q n - |_i) 

0 0 ... 0 -F(g n +i) 
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or, in a more compact form, 


^n+l X K+l (Qn -\-1 *7l) •• • ((/n+1 Qn) X |-^-n| •> 


(38) 


where we took into consideration that F(gj) = 0 if i = 1, 2,n. 
From Eq. (FT?1) one concludes that 


K+i = (<7n +i - <?i)(<?n+i - g 2 )---(^+i - q n ) ■ V n , 


so that both sides of Eq. (|38|) contain a common divisor (qq+i— <?i)(q , r)+i — q , 2 )---( ( ?n+i — <2Vi) 7^ 0. 
Consequently, A n+1 V^ = \A n \ , Q. E. D. 


A. Relationship to the inverse scattering problem 

Let us apply Theorem 1 to the alternant (1291) generated by the functions (1301) . which are 
polynomials in a variable x 2 . Using Eqs. (I34|) and (l35]h one gets 


®01 @n— l(^2j ^3 1 ••••> K-n ), 

flu @n— 2 (^2; .3■ ^n)i ••■i ®n—1,1 lj ®nl 

«02 = 0ra-l( K l> ^3; •••> K n), 

®12 @n— 2 (^1; 3■ •••) l^nji 1,2 1; ®n2 


0, 

0, 


®0n ^"n—l(^l; ^ 2) ^n—l)j 

ain ^"n— 2 (^1) 2■ l^n— l)j •••) ®n— l,n 1; ®im 0- 
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Consequently, 


D n l / n(Q , i, Q2i •••5 Qn) X |-Rn| ? l-^n 


^"n—1(^2}^3 ^"n) @n— 2(^2? ^3? • •• ? ^n) ••• 1 0 
(J n _i( aci, K3,..., ^77,) cr n —2 (fiq, K3, ..., l ^ n ) ... 1 0 

1 (^i, ^2, •..,1) 0"n—2(^1? ^2? •••) ^ n — l) ••• 1 0 

S n 5 n _x ... S! So 


^n—1(^2 5 ^3j.'.j^n) ^n—2(^2? ^3> • • • 1 ^n) ••• 1 
(J n _1 (Ki,/^3, ..., (J n _2(/"Cl, K3, ..., K n ) ... 1 

(J n _1 (/^i, K25 • • • 5 /^n—l) ®n— 2(^1? ^2? •• •? ^n— l) ••• 1 


(39) 


Here we defined a new determinant \R n \ which seems to be another alternant, so we can 
apply Theorem 1 to evaluate it. To be convinced that \R n \ is indeed an alternant, we have 
to specify the generating functions. Obviously, 


F n (x) = l = a 0 (Ki,K 2 ,...,Kn), (40) 

F n -i(x) = <Ti{ki, k 2 , ..., K n ) - X. 


Also, it is easy to prove that 


F n - 2 {x) =cr 2 - iF n _!(i) = cr 2 - a x x + x 2 , 


(41) 


where we dropped the arguments to get a more compact formula. Let us agree that if no 
arguments are explicitly given then the corresponding function depends on n arguments: 

^1 } ^2? • “ 5 ^71* 

Taking, e.g., x = we can check the validity of Eq. (ITT]) . Indeed, 

F n - 2 (ki) = a 2 - <TiKi + k\ = 
y ^ ( k i + ^2 + ••• + Kn) + Ki = 

1<71<72<71 

1(^2 5 ^3 5 •••? ^n)j 


as needed according to Eq. (139|) . Continuing in the same manner, we get the following 
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result: 


F n - 3 (x) = cr 3 - xF n _ 2 (x) = o 3 - a 2 x + a x x 2 - x 3 , 


F n -i(x ) = cr n -i - xF n —i—\ (x) = y^(-l ycrj-jX 3 -> 

3=0 


-Ffc(z) = 

3=0 


^ ^ Ctjk^F j 
3=0 


ttjk 


(—1 y^n-k-j, if j <n- k, 
0, if j > n — k. 


(42) 

(43) 


Thus | R , n | is indeed an alternant with generating functions (1341) . Consequently, according 
to Eqs. (j35]L (135j) and (PT5T) . 


Ooi 

an 

^21 • 

■ Oral 

Oo2 

«12 

a 22 . 

• O n 2 

OOra 

Oln 

®2n 

■ o nn 

^ra 

^ra-1 

Fji—2 ■ 

• ^0 


^n—1 

— O'n- 

-2 

Ora—3 

... (-ir 

-2 ( 
Ol (- 

-If - 1 

0 

®n —2 

— O n - 

-3 

O n -4 

... (-iy 

i-2 

0 

0 

3 

1 

CO 

O n - 

-4 

O n -5 

0 


0 

0 





0 


0 

0 

Ol 

-1 


0 

0 


0 

0 

1 

0 


0 

0 


0 

0 

-l) n o re ( 


Ora—l 

(-lf-Vra- 

to 

Cl 

to 


-Ol 

1 


^n—1 

~&n-2 

Ora-3 ••• ( _ 

-i \n—2 / 

-1) Oi (- 

l)ra— ! 




®n—2 

1 

Cl 

3 

CO 

On—4 ••• ( 

1 

1 

to 

0 




O n —3 

Ora—4 

Ora-5 ••• 

0 

0 







0 

0 

5 



01 

-1 

0 ... 

0 

0 




1 

0 

0 ... 

0 

0 




( 44 ) 
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which means that 



V n 


( 45 ) 


Indeed, we can transform Eq. (1441) . repeatedly using cofactor expansion in terms of the last 
column and applying the general definition 


det(A) = E <J {p) a lp 1 a 2p2--- a np n • (46) 

p 

Here the sum involves all possible permutations p = (pi,P 2 , ••• ,p n ) of the indices (1, 2, n ), 
cr(p) = (—l)^ 33 and N p is the number of pairwise interchanges needed to restore the natural 
order (1, 2, n). For example, cr(3,2,1) = —1, but a (4, 3, 2,1) = 1. As a result of the 
described operation, we obtain 

= cr(n,n- 1,..., 1) • (-l) [n/2] . 

v n 

On the other hand, a(n,n — 1,..., 1) = (—1), since N p = \nj 2], Consequently, 

\Rn\ = V n ■ (-l) [n/2] • (-l) [n / 2] = V n . 


In summary, we have obtained a very simple and universal recipe for calculating deter¬ 
minants defined by Eq. (I28]h 


Dn = (K ) 2 = n (Kj- - Ki ) 2 • 


(47) 


Combining Eqs. (l27j) and (I47|) . we can formulate a general and important result: 

Theorem 2: Let Ki, K 2 ,be arbitrary positive real numbers arranged in the ascending 
order, so that < H 2 < ■■■ < «n, and let D(k\, K 2 , K n ) be a determinant, defined by Eq. 
(ED). Then 


D(k u K2, . a ) = ]J 

l<i<j<n 


Kj — Ki 
Kj ~\~ K{ 


(48) 


where the product contains all possible combinations of the pairs (fiq, Ki) satisfying the 
condition 1 < % < j < n. 


V. SYMMETRIC REFLECTIONLESS POTENTIALS 

The excursion to theory of determinants concluded with a surprisingly simple final result. 
Indeed, as is seen from Eqs. (fT8l) . (12TT) and (I23T) . all coefficients in Eq. (1221) can be evaluated 
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with the help of Eq. (I48jh which means that the general formula for r-functions can be 
essentially simplified. For example, 


A 0 = s/D(Ki,K 2 ,..,K n ) = J ] 

l<i<j<n 


Kj — Kj 
Kj ~\~ Kj 


Ai 

1 D(Ki, K 2 , .., Ki-iKi + iK i+2 ...K n ) 

- TT 

Kj + Ki 

Ao 

\j D(K 1 ,K 2 ,..,K n ) 

11 

Kj — Ki 


(49) 


(50) 


etc. Thus Eq. (122]1 transforms to 
det (An 


t n = 


2 An 


cosh(a 0 + /do) + En 

l jjtl 


K j Ki 


Kj — Ki 


COSh(Qi + 0i) + 


(51) 


e n 

1<21<22<AT j ^2l ,22 


Kj Ki 




Kj ~\~ Kj 2 , \ 

cosh(a^ 2 + Piiiz) + ...+ 

Kj Kj 2 


e n 

<i[N/2]<N —,i[N/ 2] 


and Eq. m can be rewritten as 


Kj + 


Kj + Ki 2 


K j + K i[N/2] 

Kj Kj^ 


Kj Kj 2 


Kj — Ki^ N/2] 


COSh(d'j 1 ...j [JV/2] + ■> 


[N/2] J 


d 2 


V(x) = - 2 C—{\wt n {x)} 


(52) 


In general, as mentioned, the reflectionless potential is uniquely determined if 2 N param¬ 
eters K n and C n (n = 1, 2,..., N) are known. However, if one sets an additional constraint 


V(-x) = -V(x) 


then the number of necessary input parameters is twofold reduced. In other words, a symmet¬ 
ric reflectionless potential is uniquely determined by its N binding energies^ E n = — Ck\. 
Let us analyze this in more detail. 

Obviously, Eq. (HHf can only be symmetric if the arguments of all cosh functions are of 
the linear form ax + 6 with 6 = 0. It means, for example, that 

N 

5>*i = A> = ^2 c ij ) (53) 

2=1 1<2<J<AT 

AT 

K\ x± - k 2 x 2 - ... - k n x n = 2ki xi- f3 0 = -(3i = ^ c Xjl - j3 0 , (54) 

3 =2 
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where 


Cij = In 


Kjj + 


Kj — K{ 


(55) 


AT 


and consequently, 

2/sq Xi = Cij. 

J=2 

Using Eqs. 03 and Eq. (]23l) . we get a similar expression for any other combination KiXi. 
From Eqs. (TTUT) and (1551) we therefore obtain the following symmetricity conditions: 


exp (2 Ki Xi ) = 


C? 

2fiV 


iV 

n 


K j T Ki 


Kj — Ki 


, i = 1,2,.., N >2. 


If IV = 1 then 


c? 


exp (2 k, Xi) = — = 1. 


(56) 


(57) 


It can be easily shown that Eqs. (155]) are indeed the symmetricity conditions for r n and 
V(x). To this end, in full analogy with Eq. (1531) . we can write 

N 

K\ X\ - k 2 x 2 - ... - k n x n = ^2 Cij - /3 0 , (58) 

jV 1 

N 

-Ki Xi + K 2 X 2 - ... ~ K N X N = '^2 c 2 j ~ Po, 

0 +2 


N 

—K\ X\ + k 2 x 2 — ... + Kj\r X]y = ^ ^ c^j — Pq. 

Summing these equations, we get 

N 

-(N -2) '22K i x i = -(N-2)p 0 , 

i=1 

which coincides with Eq. (155j) . 

Remark: If iV = 2 then first two equations of the system (1551) are not linearly indepen¬ 
dent, since 

N N 

Pi — Po ~ ^2 c b = Po ~ ^2 c<2 i = /^2 = 0, 

i 2 

so that K\ X\ = k 2 x 2 . Consequently, in this (and only in this) special case Eqs. (I53|) - (l54|) 
must be treated as the actual symmetricity conditions, while Eq. (1551) still remains valid. 
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The next step is to complement Eq. (l58j) . for example, with another condition 

N N 

2 («! Xi + k 2 x 2 ) = E °ij +t: c 2j i ( 59 ) 

which is a direct conclusion from Eq. (I56H . Thus 

N 

2 (fiqxi + K2X2) - ^ xi + /c 2 o ;2 - H3X3- ~ k n x n = 

i= 1 

AT AT 

c ij+^ ^ — p 0 =/?i2, 

which means that 


cosh(«i 2 + (3 12 ) = cosh [(ki + /c 2 - «3 — ... — kjv ) 


is a symmetric function. Analogously, one can prove that any other term cosh(cq 1 j 2 ... + /di 1 i 2 ...) 
in Eq. (l5Tj) is a symmetric function as well. This in turn proves that the norming constants 
C n of a symmetric reflectionless potential are uniquely determined by the given binding 
energies E n . 


A. Few practical examples 


To illustrate the results, let us take, for example, N = 4. Then Eq. (151]) reads 


T A = 

f K 2 + fti 

\K 2 - Ki 

/c 4 + K3 

A4 - /C 3 


cosh(a 0 + Po) + ( K2 + 

\k 2 - K\ J 


/C 3 + /Ci \ f /C 4 + K\ 


K3 ~ K-l J \ K 4 — 


cosh(aq + P\)+ 


/C 3 T K -2 \ ( ^4 T ^2 


^3 — ^2) \K 4 ~ K 2 


cosh(a 3 + p 3 ) + + 

\ K 4 - K \) 
/C 3 T Kl \ ( K 4 + /Cl \ f /C 3 + K 2 


cosh(a 2 + P 2 ) + f K3 + Kl \ 
\l%3 - K] J 


K 3 + K 2 

«3 - «2 


X 


K 3 - Ki 

n 2 + /ci 
n 2 - k-i 
' k 2 T H\ 


k 4 — K\ 
K 4 + /Cl 
/C 4 — /Cl _ 
/C 3 + /Cl 


A3 - «2 y 
' /C 3 + /C 2 5 

A3 - «2 y 

/C 4 + /C 2 


/C2 — /Ci y y /C3 — /Ci y \ /C4 — /c 2 y \ /C4 — /C3 


/C4 - 4 - k 2 

A4 - «2. 

\ / /C4 + K 2 

I V /C 4 — /C 2 

\ / ^4 + 

I \/C 4 - /C 3 
/C4 + /c 3 


4 COSll(«4 + /? 4 ) + 

/c 4 - /c 3 y 
cosh(ai2 + ^12 )T 

COSh(o!i3 + ^13) + 
cosh(ai4 + /di 4 ), 


( 60 ) 
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while the symmetricity conditions, according to Eqs. (I54|h(l56]h can be given as 

' K 3 — K'2 Kb — K 2 K 4 - /t 3 ' 


/tiaq — k 2 x 2 — K3X3 — k 4 x 4 = —j 3 i = In 

-K1X1 + K2X2 - K3X3 - K4X4 = —/?2 = hi 

—Kix 1 - /c 2 a ; 2 + H3X3 - k 4 x 4 = -/ 3 3 = In 

— K1X1 — K 2 X 2 — K3X3 + K 4X3 = —04 = In 

Summing the corresponding sides of Eq. (16TT) we get 


/C 3 + K 2 K 4 4" /C2 ^4 “1“ ^3 
/C 3 — /Ci /C 4 — Ki K 4 — /C 3 
/C 3 + K\ K 4 + K 4 + K 3 
K 2 ~ K\ K4 — K\ K4 — K 2 
/C 2 + K 4 + /ti K 4 + K 2 
K 2 ~ K 1 K 3 — Ki K 3 — K 2 
K-2 + K 1 ^3 + Ki /C 3 + /t 2 


( 61 ) 


K1X1 + K 2 x 2 + K3X3 + K4X4 + A) = 0 , 


which means that 

cosh(a 0 + Po) = cosh [(/c x + k 2 + /c 3 + /c 4 ) a;] 

is a symmetric function. 

Analogously, substructing the sides of the last two equations from the corresponding sides 
of the first two equations of (16TT) . we get 

KlXi + K2X2 - K3X3 - K 4 X 4 + P \2 = 0 , 


which means that cosh(aq 2 + fti 2 ) is a symmetric function. Continuing in a similar manner, 
it is easy to be convinced that cosh(a 13 + /3 i3 ) and cosh(a 14 + p 14 ) are symmetric functions 
as well. 

Example 1: To be more specific, let xq = 1/k\ be the length unit and Eq = Ck\ , the 
energy unit (i.e., Xq — 1 and Eq = 1). The simplest and the best known symmetric potential 
then corresponds to 

K n = n, n = 1 , 2 , 3 , 4 . 


Therefore, according to Eq. (l60lh 

3 . 4.5 3-5-6 

T 4 = cosh [(ki + k 2 + k 3 + k 4 ) x] + -——— cosh [(k\ — k 2 - n 3 - k 4 ) x] + -—-—- x (62) 

_L * Zi * O _L * _L * Z, 

, r/ . , 4 • 5 ■ 7 r/ . , 5 • 6 ■ 7 

cosh [(/c 2 - K\ - k 3 - k 4 ) x\ + 211 cosh 3 - «i - K 2 - K. 4 ) x\ + ^ x 

4 . 5 . 5.6 3 . 5 . 5.7 

cosh [(/c 4 - Ki - k 2 - K 3 ) x\ + 2 - - - cosh [(fiq + k 2 - /c 3 ~ k 4 ) x\ + ^ ^ ^ x 

cosh [(«i + /c 3 — k 2 - k 4 ) x\ + ^221 COsl1 - K2 _ K s) = 

cosh (lOx) + 10 cosh ( 8 x) + 45 cosh ( 6 x) + 120 cosh (4a:) + 210 cosh ( 2 x) + 126. 
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At first sight Eq. (1621) may seem impractical. However, with the help of standard transfor¬ 
mation formulas (obtained from corresponding trigonometric formulas by replacing x —> ix ): 

cosh (lOx) = 512cosh 10 (x) — 1280cosh 8 (x) + 1120 cosh b (x) — 400cosh 4 (x) + 50cosh 2 (x) — 1, 
cosh (8x) = 128cosh 8 (x) — 256cosh 6 (x) + 160 cosli 4 (x) — 32cosh 2 (x) + 1, 
cosh (6x) = 32 cosh 6 (x) — 48 cosh 4 (x) + 18 cosh 2 (a;) — 1, 
cosh (4a;) = 8 cosh 4 (a;) — 8 cosh 2 (a;) + 1, cosh (2a;) = 2 cosh 2 (a;) — 1, 


the result is as follows: 


T 4 = 512 cosh 10 (x). 


Thus 


and 


T[ 

[lnT 4 (x)] / = —/ = 10tanh(x) 

T,i 


V(x) = 


20 


N(N+ 1) 


cosh 2 (x) cosh 2 (x) 

exactly as needed. 

Example 2: Next, let us construct a reflectionless approximant to a symmetric rectangular 
potential with four energy levels (see Fig. [[]). These levels can be determined from (see, 
e.g.,—, Sect. 11.9) 


Ki Xi tan (m Xi) = -(^Xj) 2 , 


C 


Kji Xi 


U 0 a 2 


(^2 %i) t 


(63) 

(64) 


tan (Ki x^ V C 

where a and Uq denote the half-width and the depth of the potential well, respectively. 
Eq. (163|) Exes the symmetric and (1641) - the antisymmetric solutions to the Schrodinger 
equation. Again, it is convenient to use dimensionless units for the length and energy, 
taking x 0 = a = 1 and E 0 = C/a 2 = 1. In addition, let us fix 


U 0 a 2 

~C~ 


= 5. 


Then the system has four discrete levels (as assumed) corresponding to 


K\ = 1.3064400089, 

k 2 = 2.5957390789, (65) 

= 3.8374671080, 
k 4 = 4.9062951521. 
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Dimensionless coordinate 


FIG. 1. The isospectral reflectionless approximant (red dashed line) to a symmetric rectangular 
potential (blue solid line). The common energy levels are shown by green dashed lines. 


In this case Eq. (I60p cannot be further simplified, but this is not any serious problem. 
Indeed, let us define the coefficients and B { , such that 


T 4 = Ai cosh (Bix). 

i 


Then the corresponding potential becomes 
V{x) = -2 C 


f Y A t B‘f cosh (Bi x) 

A,B, sinh (Bix)~ 

i 

2 1 

j ^2 Ai cosh (Bi x) 

Y^Ai cosh {Bi x) 

i 

i 


( 66 ) 


The result for the input data (IS5j) can be seen in Fig. |T| 

Example 3: Fig. [2] demonstrates two isospectral potentials corresponding to the following 
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FIG. 2. The isospectral reflectionless substitute to a Morse potential (see the explanations to 
Example 3). The common energy levels are shown by green dashed lines. 


set of input parameters: 


Ki = 1/8, k 2 = 3/8, k 3 = 5/8, k 4 = 7/8, 


(67) 


the four energy levels being E n = —Ck\ (n = 1,2, 3,4) as previously. The solid curve in 
this figure corresponds to a Morse potential 1 


V{x) 


= exp 


2 ax 

x 0 


, ax\ 

- 2exp '--) 


( 68 ) 


taking xq = \Jh 2 /(2mD) = 1, E 0 = D = 1, and consequently, C = 1. The energy eigenvalues 
read 

E n = -D{ 1 - W+1/2 ) 2 , (69) 

a 

where a = y/D/C/a. In Fig. [2l the value a = 4 has been chosen, so that a = 1/4 in our 
dimensionless units. As in the previous example, the red dashed curve shows the symmetric 
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reflectionless potential derived by Eq. (I66|) from the input parameters (1671) . 


VI. CONCLUSION 

The main result of this work is a general formula for calculating r-functions. This im¬ 
portant formula, Eq. (l5Tjh is a direct conclusion from Theorem 2 that has been proved 
with the help of well-known methods of the theory of determinants. We demonstrated that 
r-functions can be expanded in terms of special determinants called alternants^. Any alter¬ 
nant has a divisor - the Vandermonde’s determinant of the same order, while the quotient 
can be uniquely expressed as a polynomial in elementary symmetric functions (fool) (see The¬ 
orem 1). Moreover, in the case of alternants related to the inverse scattering problem this 
quotient equals unity, i.e., the alternant itself equals the Vandermonde’s determinant. These 
useful properties of alternants is the key to a very simple final result expressed by Eq. (ED- 

Using Eqs. (l5Tf - (l52]h one can uniquely reconstruct any reflectionless one-dimensional 
potential on the full line (—oo < x < oo), provided that the 2 N input parameters K n and C n 
(n = 1,2, ...,N) are known. Moreover, if the result is expected to be a symmetric function 
of the coordinate x then the problem can be uniquely solved on the basis of the N binding 
energies E n = —Ck^. Compared with the direct use of Eq. (J9]), the described approach 
significantly reduces computational efforts. Indeed, the expansions (122]) and ()5T]) contain 
only 2 n ~ 1 members, while Eq. d^J) requires evaluation of a determinant with N\ members. 

The efficiency of the method has been explicitly demonstrated for the case N = 4, and 
there is no doubt that the algorithm can be successfully applied to a much higher number 
(in principle, to an arbitrary number) of given binding energies. The described approach 
can also be applied to building TV-soliton solutions to the Korteweg-de Vries equation, but 
this would be a subject for another paper. 
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